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Abstract 

Let X be a set of points whose coordinates are known with limited 
accuracy; our aim is to give a characterization of the vanishing ideal 2T(X) 
independent of the data uncertainty. We present a method to compute a 
polynomial basis B of X(X) which exhibits structural stability, that is, if X 
is any set of points differing only slightly from X, there exists a polynomial 
set B structurally similar to B, which is a basis of the perturbed ideal X(X). 
Keywords: Empirical points, vanishing ideal, border bases. 
MSC: 13P10, 65F20, 65G99. 

1 Introduction 

In this paper we present a method for computing "structurally stable" border 
bases of ideals of points whose coordinates are affected by errors. 

If X is a set of "empirical" points, representing real-world measurements, 
then typically the coordinates are known only imprecisely. Roughly speaking, 
if X is another set of points, each differing by less than the uncertainty from 
the corresponding element of X, then the two sets can be considered as equiva- 
lent. Nevertheless, it can happen that their vanishing ideals have very different 
bases — this is a well known phenomenon in Grobner basis theory. In order to 
emphasize the "numerical equivalence" of X and its perturbation X, we look for 
a common characterization of the vanishing ideals I(X) and X(X). More pre- 
cisely our goal is to determine a polynomial basis B of the vanishing ideal 2~(X) 
which exhibits structural stability: namely, there is a basis B for the perturbed 
ideal I(X), sharing the same structure as B, and whose coefficients differ only 
slightly, provided that X differs from X by only a small amount (up to some 
limit). 

The decision to use border bases to describe vanishing ideals of sets of em- 
pirical points was due to two main reasons: border bases have always been 
considered a numerically stable tool (see [9], HB); furthermore, it is easy to 
study their structure, i.e. the support of their polynomials, as it is completely 
determined once a suitable order ideal O has been chosen. 
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We introduce the notion of stable quotient basis: given a set X of empirical 
points and a permitted tolerance e, a stable quotient basis O guarantees the 
existence of an O-border basis B for the vanishing ideal I(X) where X is any 
set of points perturbed by amounts less than the tolerance e. Once a stable 
quotient basis O has been found, the corresponding stable border basis can be 
obtained by some simple combinatorical and linear algebra computations; so we 
focus our attention on determining O. 

An alternative approach to the problem, presented in [7], is to use singular 
value decomposition of matrices to obtain a set of polynomials which are not 
required to vanish on X but must nevertheless assume particularly small values 
there. In contrast, a stable border basis always comprises polynomials which 
vanish on X. 

This paper is organized as follows. In Section 2 we introduce the concepts 
and tools we shall use. Section 3 provides a formal description of our prob- 
lem. The main result, the SOI algorithm for computing a stable order ideal, 
is presented in Section 4. In Section 5 we give some numerical examples il- 
lustrating the functioning of the algorithm. Finally, Section 6 is an Appendix 
which contains the proof of a basic result about the first order approximation of 
rational functions, useful for the error analysis of the sensitivity of the border 
basis computation. 

2 Basic definitions and notation 

This section contains basic definitions and notation used later in the paper. 
To simplify the presentation, we shall implicitly suppose that each finite set of 
points or polynomials is in fact a tuple, so that the elements are ordered in some 
way, and we can refer to the fc-th element using the index k. 

Let n > 1; we recall (see [HUH]) some basic concepts related to the polynomial 
ring P = R[xi, . . . ,»„]. 

Definition 2.1 Let X = {pi, . . . ,p s } be a non-empty finite set of points ofW 1 
and let G = {gi, . . . , be o, non-empty finite set of polynomials. 

(a) The ideal I(X) = {/ € P \ fipi) = Vp,; € X} is called the vanishing 
ideal o/X. 

(b) The M.-linear map evalx : P — > R s defined by evalx(J) = (f(pi ),-■-, f(Ps)) 
is called the evaluation map associated to X. For brevity, we write /(X) 
to mean evalx (/)• 

(c) The evaluation matrix ofG associated to X, written Mq(X) S Mat* fc(R), 

is defined as having entry equal to gj(pi), i.e. whose columns are the 
images of the polynomials gj under the evaluation map. 

Definition 2.2 Let T™ be the monoid of power products of P and let O be a 

non-empty subset o/T™. 
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(a) The factor closure (abbr. closure,) of O is the set O of all power products 
in T" which divide some power product of O. 

(b) The set O is called an order ideal if O = O, i.e. if O is factor closed. 

(c) Let I C P be a zero- dimensional ideal, and s = dim(P/J); if O is factor 

closed and the residue classes of its elements form a basis of Pj I then we 
call it a quotient basis for I. 

(d) Let O be an order ideal; the border dO of O is defined by 

dO = ( Xl O U ... U x n O) \ O 

(e) If O is an order ideal then the elements of the minimal set of generators of 

the monomial ideal corresponding to T n \0 are called the corners of O. 

Definition 2.3 Let O = {tx, . . . , t^} be an order ideal, and let dO = {bi, . . . , b„} 
be the border of O. Let B = {g\, . . . , g u } be a set of polynomials having the form 
9j = — Y]j—i ctijti where each aij € K. Let I C P be an ideal containing B. 
If the residue classes of the elements of O form a H-vector space basis of P/I 
then B is called a border basis of I founded on O, or more briefly B is an 
O-border basis of I. 

Proposition 2.4 (Existence and Uniqueness of Border Bases) 

Let I C P be a zero-dimensional ideal, and let O = {t\, . . . ,t^} be a quotient 
basis for I. Then there exists a unique O-border basis B of I. 

Proof: See Proposition 6.4.17 in [S]. 

Later on, in order to measure the distances between points of R", we will 
use the euclidean norm || • ||. Additionally, given an n x n positive diagonal 
matrix E, we shall also use the weighted 2-norm || • \\e as defined in [5]. For 
completeness, we recall here their definitions: 



and \\v\\ E :=\\Ev\\ 

We recall the definition of empirical point (see [TT], [5]). 

Definition 2.5 Let p € R n be a point and let e = (ei, . . . ,s n ), with each Si € 
M + . be the vector of the componentwise tolerances. An empirical point p e is 
the pair (p,e), where we call p the specified value and e the tolerance. 

Let p 6 be an empirical point. We define its ellipsoid of perturbations: 

N(p £ ) = {peR n : \\p-p\\ E <l} 

where the positive diagonal matrix E — diag(l/e , i, . . . , l/e n ). This set con- 
tains all the admissible perturbations of the specified value p, i.e. all points 
differing from p by less than the tolerance. 
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We shall assume that all the empirical points share the same tolerance e, 
as is reasonable if they derive from real-world data measured with the same 
accuracy. In particular this assumption allows us to use the E'-weighted norm 
on 1" to measure the distance between empirical points. 

Given a finite set X £ of empirical points all sharing the same tolerance e, we 
introduce the concept of a slightly perturbed set of points X by means of the 
following definition. 

Definition 2.6 Let X e = {jp\ , . . . ,p e s } be a set of empirical points with uniform 
tolerance e and with X C K™. Each set of points X = {p{, . . . ,p s } C M" whose 
elements satisfy 

s 

(pi,...,p s )el[N(pf) 

i=l 

is called an admissible perturbation o/X e . 

Finally we introduce the definition of distinct empirical points. 

Definition 2.7 The empirical points p\ and p\, with specified values pi,P2 € 
W l , are said to be distinct if 

N( P l) n N(pl) = 
3 The formal problem 

We shall use the concept of empirical point to describe formally the given un- 
certain data: the input X is viewed as the set of specified values of X £ , which 
consists of s distinct empirical points all sharing the same fixed tolerance e. 

Given the set X, we want to determine a numerically stable basis B of the 
vanishing ideal I(X). Intuitively, a basis B of X(X) is considered to be struc- 
turally stable if, for each admissible perturbation X of X £ , it is possible to 
produce a basis B of I(X) only by means of a slight and continuous variation of 
the coefficients of the polynomials of B, that is if there exists a basis B of I(X) 
whose polynomials have the same support as the corresponding polynomials 
of B. Given a polynomial basis B, we will call the union of the supports of its 
polynomials the structure of B. 

A good starting point for us is the concept of border basis (see [5], [TT]). 
In fact the structure of a border basis is easily computable and completely 
determined by the quotient basis O upon which the border basis is founded (see 
Definition 12. 3p . Using border bases, the problem of computing a structurally 
stable representation of the vanishing ideal I(X) thus reduces to the problem of 
finding a quotient basis O for X(X) valid for every admissible perturbation X. 
The following deinition captures this notion and generalizes it to any order ideal. 

Definition 3.1 Let O be an order ideal, then O is stable w.r.t. X e if the eval- 
uation matrix Ma (X) has full rank for each admissible perturbation X of X e . 
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The following proposition highlights the importance of stable quotient bases. 

Proposition 3.2 Let X e be a set of s distinct empirical points, and let O = 
{t\,...,t s } be a quotient basis for X(X) which is stable w.r.t. X s . Then, for 
each admissible perturbation X of X s , the vanishing ideal X(X) has an O -border 
basis. Furthermore, if dO = {b\, . . . , b v } is the border of O then B consists of v 
polynomials of the form 

s 

9o = b 3 <~ ^2 Vijti for j = 1 ... 1/ (1) 

i=l 

where the coefficients £ R satisfy 

s 

& j (X) = ^o ij -t i (X) 
i=i 

Proof: Let X be an admissible perturbation of X £ and let evalj : P — > R* 
be the R-linear evaluation map associated to the set X. It is easy to prove 
that I(X) = ker(eval^) and consequently, that the quotient ring P/X(X) is 
isomorphic to R s as a vector space. Since O is stable w.r.t. the empirical 
set X e , it follows that {ti(X), . . . , t s (X)} are linearly independent vectors. More- 
over #X = #0, so the residue classes of the elements of O form a R- vector space 
basis of P/X(X) L 

Let vj = bj (X) be the evaluation vector associated to the power product bj 
lying in the border dO; each vj can be expressed as 

s 

Vj = ^ dijti(X) for some G R 

For each j we define the polynomial gj — bj — J2t=i a ijU\ by construction 
eval^(gj) = 0, and so B = {gi, . . . , g v } is contained in X(X); it follows that B 
is the O-border basis of the ideal X(X) . 

We observe that the coefficients a^- of each polynomial gj £ B are just the 
components of the solution ctj of the linear system Mo(X) ■ ctj — bj(X). It 
follows that ctij are continuous functions of the points of the set X and so, since 
the order ideal O is stable w.r.t. X e , they undergo only continuous variations as 
X changes. Now, the definition of stable border basis follows naturally. 

Definition 3.3 Let X s be a finite set of distinct empirical points, let O be a 
quotient basis for the vanishing ideal X(X). Lf O is stable w.r.t. X e then the 
O-border basis B forX(X) is said to be stable w.r.t. the set X e . 

The problem of computing a stable border basis of the vanishing ideal of a 
set X £ of empirical points is therefore completely solved once we have found a 
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quotient basis O which is stable w.r.t. X £ . If O exhibits these characteristics, 
Proposition 13.21 and the subsequent observation on the continuity of the coef- 
ficients ctij prove the existence of the corresponding stable border basis of the 
ideal X(X) . The problem of the effective computation of a stable quotient basis 
is addressed in section 0] 

We end this section by observing that any O-border basis of the vanishing 
ideal Z(X) is stable w.r.t. X s for a sufficiently small value of the tolerance S. 
This is equivalent to saying that any quotient basis O of Z(X) has a "region of 
stability" , as the following proposition shows. 

Proposition 3.4 Lei X be a finite set of points o/K™ andl(X) be its vanishing 
ideal; let O be a quotient basis for 2(H). Then there exists a tolerance S — 
(6i, ■ ■ ■ , 5 n ), with Si > 0, such that O is stable w.r.t. X 5 . 

Proof: Let Mo(X) be the evaluation matrix of O associated to the set X; 
then Mq(X) is a structured matrix whose coefficients depend continuously on 
the points in X. Since, by hypothesis, the O-border basis of the vanishing 
ideal I(X) exists, it follows that Mo(X) is invertible. Recalling that the deter- 
minant is a polynomial function in the matrix entries, and noting that the entries 
of Mo(X) are polynomials in the points' coordinates, we can conclude that there 
exists a tolerance 5 = (<5i, . . . , S n ), with each Si > 0, such that det(Mo(X)) ^ 
for any perturbation X of X. 

Nevertheless, since the tolerance e of the empirical points in X e is given a 
priori by the measurements, Proposition 13.41 does not solve our problem. If the 
given tolerance e is larger than the "region of stability" of a chosen quotient 
basis O, the corresponding border basis will not be stable w.r.t. X e ; such a 
situation is shown in the following example. 

Example. Let X e be the set of empirical points having 

X={(-l,-5), (0,-2), (1,1), (2,4.1)} CM 2 

as the set of specified values and e = (0.15, 0.15) as the tolerance; let 

X={(-l + ei,-5 + e 2 ), (e 3 ,-2 + e 4 ), (l + e 5 ,l + e 6 ), (2 + e 7 , 4.1 + e 8 )} 

be a generic admissible perturbation of X e , where the parameters € K satisfy 

||(ei,e 2 )|U < 1 |Ke 3 ,e 4 )|| E <l || (e 5 , e 6 )|| B < 1 ||(e 7 , e 8 )|| E < 1 

Consider first 0\ = {1, y, x, y 2 }, which is a quotient basis for I(X). The corre- 
sponding border basis B\ of X(X) is not stable w.r.t. X e . Indeed, consider the 
perturbation X = {(-1,-5), (0,-2), (1,1), (2,4)}j)fX £ . The evaluation ma- 
trix Mfjj (X) is singular, so no Oi-border basis of X(X) exists. It follows that 0\ 
is not stable w.r.t. X e since its "region of stability" is too small w.r.t. the given 
tolerance e. 

Now consider the quotient basis 02 = {1, y, y 2 , y 3 }, which is stable w.r.t. X e . 
In fact, for each perturbation X of X e , we see that Mo 2 (X) is a Vandermonde 
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matrix whose determinant is equal to (e 4 — e2 +3)(e6 — &2 + 6)(es — &2 +9.1) (eg — 
e 4 + 3)(e 8 — e 4 + 6.1)(e 8 — e 6 + 3.1). Since each |ej| < 0.15, it follows that, for each 
perturbation X, the matrix Afe> 2 (X) is invertible, and so it is always possible to 
compute an 02-border basis of the ideal Z(X) . <J) 



4 A practical solution 

In this section we address the problem of computing an order ideal O stable 
w.r.t. X e , a finite set of distinct empirical points, and also the corresponding 
stable border basis when it exists. 



Since in real-world measurements the tolerance e present in the data is rel- 
atively small, our interest is focused on small perturbations X of the empirical 
set X e . For this reason our approach is based on a first order error analysis of 
the problem. We present in Section 14.31 an algorithm which computes a stable 
order ideal O (up to first order). In order to investigate the stability of the 
order ideal O we use some results on the first order approximation of rational 
functions (see Section 14. ip and we introduce a parametric description of the 
admissible perturbations X of X e (see Section 14.21) . 

If the output of the algorithm is actually a quotient basis then the cor- 
responding stable border basis B exists for I(X). To determine B it suffices 
to find the border of O (a simple combinatorical computation), and then for 
each element of the border solve a linear system (as described in the proof of 
Proposition 13. 2p . 



4.1 Remarks on first order approximation 

Let n £ N; let F = K(ei, . . . , e n ) be the field of rational functions on K and 
let / £ F. We use multi-index notation to give the Taylor expansion of / in a 
neighbourhood of 



/(d, . . . ,e„) = ^ 



D a f(0) ca 



We recall that given a = (ai, . . . , a n ) £ N n , we have \a\ — ol\ + . . . + a n 
and a\ = ai \...a n \ Similarly D a = D* 1 ...£>£» (where = d 3 /del) and 
e — c 1 . . . t n . 

Each f £ F can be decomposed into components of homogeneous degree in the 
following way: 

. v-, , , D a f(0) a 

f = 2^ fk where /* = 2^ ~^^ e 

fc>0 \a\=k 

and where, by convention, £)( - )/ = /. Each polynomial fk £ R[ei, . . . , e n ] is 
called the homogeneous component of / of degree k. 

Analogously, we can decompose a matrix M £ M&t riS (F) into "homoge- 
neous" parts in the following way. 
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Definition 4.1 Let M = (rriij) be a matrix in Mat riS (F); the matrix = ((rriij)k), 

where (rriij)k = J2\ a \=k ~ "j^^ e" e K[ei, . . . , e„], is called the homogeneous 
component of M of degree k. 

The following proposition characterizes the homogeneous components of de- 
grees and 1 of the solution and residual of a least squares problem. 

Proposition 4.2 Let r, s e N with r > s; let M be a matrix in M&t r ^ s (F) 
and let v be a vector in Mat r! i(F). Let x € Mat Si i(.F) and p € M&t r ,i(F) be 
respectively the solution and the residual of the least squares problem Mx ~ v. 
The homogeneous components of degrees and 1 of x are 

x = {MlM^Mlvo (9] 
x x = (MqMq)~ 1 (MqVi + M[v - MqMiXq — M{M x Q ), [Z) 

Moreover, the homogeneous components of degrees and 1 of p are 

pa = v - M xq .g. 
p 1 = ui - M x 1 - Mix 1 ' 

Proof: This lengthy proof has been deferred to an appendix. 

Since the residual p is orthogonal to the columns of the matrix M, we have 

MqP = and M\po + M^pi = 

but this does not imply that the vector p + p\ is orthogonal to the columns of 
M + Mi. Nevertheless, since 

(M + Mi)*(p + Pi) = M*p + M* Pl + Mlpo + M{p x = M{p x 

we can assert that the vector p + p\ is orthogonal to the columns of M + M\ , 
up to first order. 

4.2 A parametric description of X £ 

Let X e = {p\ , . . . ,p £ s } be a finite set of distinct empirical points with specified 
values X C 1"; we represent an admissible perturbation X = {p\, . . . ,p s } of X e 
by using first order infinitesimals for the perturbation in each coordinate. In 
particular we express X as a function of ns variables 

e = (en , . . . , e s i, ei2, ■ ■ ■ , e S 2, ■ ■ ■ , e\ n , . . . , e sn ) 

called error variables; specifically, we have 

Pk = (Pki + efei, p k 2 + ek2, ■ ■ -Pkn + e kn ) 

The condition that each p^ is an admissible perturbation of the point p^ is 
equivalent to the following: 

||(efci,... ) e fcn )|| B <l (4) 
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We denote by X(e) = (pi(e), . . . ,p s (e)) a generic admissible perturbation 
of X e . We observe that the coordinates of each perturbed point Pfc(e) are el- 
ements of the polynomial ring R — M[e] and that each variable ejy represents 
the perturbation in the j-th coordinate of the original point p/.. The domain of 
the perturbed set X(e), viewed as a function of ns variables, is denoted by D e . 
Obviously, if e e D e we have 

n s n 
3=1 k=l ] = 1 

and consequently 

l|e||<V5|N (5) 

To keep evident the dependence on the error variables e, we extend the 
concepts of evaluation map of a polynomial / G P and evaluation matrix of 
a set of polynomials G = {g%, . . . , (?/-} C P (see Definition 12. ip to a generic 
perturbed set X(e), using the following notation: 

eval s (e)(/) = (/(pi(e)), . . . , /(p a (e))) eRx...xR = R s 

for brevity denoted by /(X(e)); similarly we write the evaluation matrix 

M G (X(e)) = (. 9l (X(e)),...,.g fc (X(e))) 

4.3 The SOI Algorithm 

In this section we present the SOI algorithm which computes, up to first order, 
an order ideal O stable w.r.t the empirical set X e . 

The strategy for computing a stable order ideal O is the following. As in 
the Buchberger-Moller algorithm ([3], pQ) the order ideal O is built stepwise: 
initially O comprises just the power product 1; then at each iteration, a new 
power product t is considered. If the evaluation matrix Af OU { t -j(X(e)) has full 
rank for all e in D £ then t is added to O; otherwise t is added to the corner set 
of the order ideal. 

A first observation concerns the choice of the power product t to analyze 
at each iteration: any strategy that chooses a term t such that the set O U 
{t} preserves the property of being an order ideal can be applied. A possible 
technique is the one used in the Buchberger-Moller algorithm, where the power 
product t is chosen according to a fixed term ordering a. The version of the 
SOI Algorithm presented below employs this latter stategy. Note that a is only 
used as a computational tool for choosing t; in fact the final computed set O is 
not, in general, the same as that which would be obtained processing the set X 
by the Buchberger-Moller algorithm with the same term ordering. 

Another observation concerns the main check of the algorithm: note that 
the rank condition is equivalent to checking whether p(e), the component of 
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the evaluation vector f(X(e)) orthogonal to the column space of the matrix 
Mo(X(e)), vanishes for any e G D e . This check is greatly simplified by our 
restriction to first order error terms. 

Algorithm 4.3 (Stable Order Ideal Algorithm) 

Let a be a term ordering on T" and let X e = {p\ , . . . ,p e s } be a finite set of 
distinct empirical points, with X C M n and a common tolerance e — {e\, . . . , £„). 
Let e = (en, . . . , e sn ) be the error variables whose constraints are given in Q. 
Consider the following sequence of instructions. 

51 Start with the lists O = [1], L = [x\, . . . , x n ], the empty list C = [ ], 

and the matrices Mq £ Mat Sj i(M) with all the elements equal to 1, and 
Mi € Mat S; i(i?) with all the elements equal to 0. 

52 If L — [] then return the set O and stop. Otherwise let t = min cr (_L) and 

delete it from L. 

53 Let vq and v\ be the homogeneous components of degrees and 1 of the 

evaluation vector v = t(X(e)). Solve up to first order the least squares 
problem Mo(X(e)) a ( e ) ~ v > by computing the vectors 



54 Let Ct € Mat S)Sn (R) be such that p\ — Cte. Compute the minimal 2-norm 

solution e of the underdetermined system Cte = — po JEj- 

55 If ||e|| > -\/s|l e ll then adjoin the vector vq as a new column of Mo and the 

vector vi as a new column of Mi. Append the power product t to O, and 
add to L those elements of {x\t, . . . , x n t} which are not multiples of an 
element of L or C. Continue with step S2. 

56 Otherwise append t to the list C , and remove from L all multiples of t. 

Continue with step S2. 

Theorem 4.4 Algorithm \4-3\ stops after finitely many steps and returns a set 
O C T" which is an order ideal stable (up to first order) w.r.t. the empirical 
set X s . Furthermore, if #0 = s then X(X) has a corresponding stable border 
basis w.r.t. X e . 

Proof: First we claim that the vectors po, p\, ot\ computed in step 53 
are the homogeneous components of degrees and 1 of the residual p(e) and of 
the solution a(e) to the least squares problem 



po = v - M a 
pi = Vi- M ai 



Mia 



where 



a = (M^Mo^M'vo 

a t = (MqM )^ 1 (MqV 1 + M{vo 



Mi 



lM ia o - M{Moao). 



M (X(e)) a(e) » i(X(e)) 



(6) 
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where t is the power product being considered at the current iteration, and O 
is the order ideal computed so far. To prove this claim it is sufficient to apply 
Proposition 14.21 to ^ and to observe that the matrices Mo and Mi coincide 
with the homogeneous components of degrees and 1 of Mo(X(e)). Clearly, 
this is true at the first iteration, since Mgi(X(e)) = (1 . . . 1)*. We apply induc- 
tion on the number of iterations. Assume that Mq and M\ are the components 
of degrees and 1 of Mo(JL(e)) and suppose that the power product r is added 
to O. Since the last column of Mo U { r }(X(e)) is given by r(X(e)), whose com- 
ponents of degrees and 1 are r$ and r\ respectively, the new matrices [Mo, To] 
and [Mi, r{\ are the components of degrees and 1 of M OU { r y(X(e)). We con- 
clude that the vectors po + pi(e) and ao + ai(e) coincide with p(e) and a(e), 
up to first order. 

Now we prove the finiteness and the correctness of Algorithm 14.31 
First we show finiteness. At each iteration the algorithm performs either step S5 
or step S6. We observe that step S5 can be executed at most s — 1 times; in 
fact, when M becomes a square matrix, i.e. after s — 1 iterations of step S5, 
the residual vector po is zero, and consequently the minimal 2-norm solution e 
of the linear system Cte = —po is also zero. Moreover, step S5 is the only place 
where the set L is enlarged with a finite number of terms, while each iteration 
removes from L at least one element; we conclude that the algorithm reaches 
the condition L = [ ] after finitely many iterations. 

In order to show correctness we prove, by induction on the number of iter- 
ations and using a first order error analysis, that the output set O is an order 
ideal stable w.r.t. X e . This is clearly true after zero iterations, i.e. after step SI 
has been executed. By induction assume that a number of iterations has al- 
ready been performed and that the set O satisfies the given requirements; let us 
follow the steps of the new iteration, in which a power product t is considered. 
If step S6 is performed the claim is true because O does not change. Otherwise, 
if step S5 is performed, the set O* — O U {t} is an order ideal by construction. 
Further, since the minimal 2-norm solution e of the linear system Cte = —po 
satisfies condition ||e|| > yfs \\e\\ it follows that e does not belong to D e and that 
the vector po + pi(e) does not vanish as e varies in D e . Therefore, up to first or- 
der, we can consider p(e) as a non- vanishing vector for each perturbation X(e), 
i.e. the matrix Me>*(X(e)) has full rank for each e E D e . 

For the last part of the theorem we simply observe that when #0 = s then 
O is a quotient basis; the rest is immediate. 

5 Numerical examples 

In this section we present some numerical examples to show the effectiveness 
of the SOI algorithm. Our algorithm is implemented using the CH — V language 
and the CoCoALib, see [3], and all computations have been performed on an 
Intel Pentium M735 processor (at 1.7 GHz) running GNU/Linux. In all the 
examples, the SOI algorithm is performed using a fixed precision of 1024 bits 
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for the RingTwinFloat implemented in CoCoALib, and the degree lexicographic 
term ordering a; in addition, the coefficients of the polynomials are displayed 
as truncated decimals. 

The first two examples show how the SOI algorithm detects the simplest 
geometrical configuration almost satisfied by the empirical set X e . 

Example. Almost aligned points 

We consider the empirical set X e given in Example we recall here the points 
in X 

X={(-l,-5), (0,-2), (1,1), (2,4.1)} CM 2 
and the tolerance e = (0.15, 0.15). 

Applying algorithm SOI to X e we obtain the quotient basis O = {1, y, y 2 , y 3 } 
which is stable w.r.t. X e , as we proved in Example [31 As O is a quotient basis 
we can compute the border basis founded on it: 

x + 0.0002?/ 3 + 0.0012?/ 2 - 0.3328?/ - 0.6686 

xy + 0.0008?/ 3 - 0.3286?/ 2 - 0.6643?/ - 0.0079 

B { xy 2 - 0.3301?/ 3 - 0.6471?/ 2 + 0.0098?/ - 0.0326 

xy 3 - 0.0199?/ 3 - 7.1199?/ 2 - 7.3933?/ +13.533 

y 4 + 1.9y 3 - 21. 6?/ 2 - 22.3?/ + 41 

Note that the lowest degree polynomial of B, f = x + 0.0002?/ + 0.0012?/ 2 - 
0.3328?/ — 0.6686, highlights the fact that X contains "almost aligned" points. 
In fact, if we neglect the terms with smallest coefficients, / simplifies to a; — 
0.3328?/ — 0.6686. Since the coefficients of a polynomial are continuous functions 
of its zeros and the quotient basis O is stable w.r.t. X e , we can conclude that 
there exists a small perturbation X of X containing aligned points and for which 
the associated evaluation matrix Mq(X) is invertible. A simple example of such 
a set is given by X= {(-1,-5), (0,-2), (1,1), (2,4)}. 

A completely different result is obtained by applying to the set X the Buchberger- 
Moller algorithm w.r.t. the same term ordering a. The c-Grobner basis Q 
of 1(X) is: 

x 2 - l/9y 2 - 121/30a; + 9/10y+ 101/45 
xy - 1/3?/ 2 -41/10X + 7/10?/ + 41/15 
y 3 + 6?/ 2 + 516243/ 100x- 171781/100?/- 172581/50 

and the associated quotient basis is O a (l(X)) = T 2 \LT CT {X(X)} = {l,y,x, y 2 }. 
We observe that O rT (l(X)) is not stable (see Example EJ) because the evaluation 
matrix Mq o (X) is singular for some admissible perturbations of X. In particular, 
the information that the points of X are "almost aligned" is not at all evident 
from Q . 

Example. Empirical points close to an ellipse 

Let X C I 2 be a set of points created by perturbing by less than 0.1 the 
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coordinates of 10 points lying on the ellipse x 2 + 0.25y 2 — 25 = 0, 

X = {(-5.07, 0.02), (4.98, 0), (3.05, 8.07), (3.01, -8.02), (-3.02, 7.99), 
(-2.98, -8), (4.01, 5.94), (3.98, -6.06), (-3.92, 6.03), (-4.01, -6)} 

Let X 6 be the set of empirical points whose set of specified values is X and whose 
common tolerance is e = (0.1, 0.1). 

Applying SOI on X e we obtain, after 11 iterations, the stable quotient basis 

= {l,y,x, y 2 , xy, y 3 , xy 2 , y A , xy 3 , xy 4 } 

We use linear algebra to compute the corresponding stable border basis B 
of X(X) . We can identify the "almost elliptic" configuration of the points of X 
by looking at / the lowest degree polynomial contained in B: 

f = x 2 + 0.273y 2 - 25.250 + 10~ 2 (0.004xy 4 + 0.020a;y 3 - 0.034?/ 4 - 0.489:n/ 2 
-0.177y 3 - 1.371xy + 9.035x + 9.810?/) 

We observe that / highlights the fact that X contains points close to an ellipse. 
In fact, if we neglect the terms with smallest coefficients, / simplifies to x 2 + 
0.273y 2 — 25.250. Since the coefficients of a polynomial are continuous functions 
of its zeros and the quotient basis O is stable w.r.t. X e , we can conclude that 
there exists a small perturbation X of X containing points lying on an ellipse 
and such that the associated evaluation matrix Me>(X) is invertible. A simple 
example of such a set is given by 

X = {(-5,0), (5,0), (3, 8), (3, -8), (-3, 8), 

(-3, -8), (4, 6), (4, -6), (-4, 6), (-4, -6)} 

Example. Empirical points close to a circle 

In this example we show the behaviour of the SOI algorithm when applied to 
several sets of points with similar geometrical configuration but with different 
cardinality. 

Let Xi,X2,X3,X4 C K 2 be sets of points created by perturbing by less than 
0.01 the coordinates of 8, 16, 32 and 64 points lying on the circumference x 2 + 
y 2 — 1 = 0, and let e = (0.01, 0.01) be the tolerance. We summarize in Table Q] 
the numerical tests performed by applying the SOI algorithm to the empirical 
set X?, for i = 1 . . . 4. The first two columns of the table contain the name 
of the processed set and the value of its cardinality. The column labelled with 
"Corners" refers to the set of corners of the stable order ideal computed by the 
algorithm; the column labelled with "Time" contains the time taken to compute 
the quotient bases. 

Note that the set of corners of the stable quotient bases computed by the 
SOI algorithm always contain the power product x 2 : this means that there is 
an "almost linear dependence" among the power products {1, y, x, y 2 , xy, x 2 } 
and that some useful information on the geometrical configuration of the points 
could be found. 
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Input 


#x« 


Corners 


Time 


Xi 


8 


{a; 2 , xy 3 ,y b } 


0.5 s 


x 2 


16 


{x 2 ,xy 7 ,y 9 } 


8.5 s 


x 3 


32 


{x 2 ,xy 15 ,y 17 } 


79 s 


x 4 


64 


{x 2 ,xy 31 ,y 33 } 


2320 s 



Table 1: SOI on sets of points close to a circle 



Example. Empirical points close to an hyperbola, a circle and a cubic 

The numerical tests suggest that in most cases the SOI algorithm computes 
a stable quotient basis, allowing us to determine a stable border basis of I(X). 
Nevertheless, this is not true in general, as the following example illustrates. 
Let X e be the set of distinct empirical points having 

X = {(1, 6), (2, 3), (2.449, 2.449), (3, 2), (6, 1)} C M 2 

as the set of specified values and e — (0.25, 0.25) as the tolerance. 
Applying the algorithm SOI to X £ , we obtain the stable order ideal O = 
{1, y, x, y 2 }; however, this is not a quotient basis, so we cannot obtain a corre- 
sponding stable border basis. This is due to the fact that the points of X lie 
close to the hyperbola xy — 6 = 0, the circle (x — 6) 2 + (y — 6) 2 — 25 = and 
the cubic y 3 — 12y 2 + 6x + 47 y — 73 = 0. So, if the tolerance e is too big, they 
"almost satisfy" all of them. 

Observe how the problem does not arise if we use a smaller tolerance, e.g. 
5 = (0.2,0.2). Applying SOI to X* 5 we obtain the stable quotient basis O' = 
{1, y, x, y 2 , y 3 }, and its corresponding border basis: 



B' = 



xy + 0.0047?/ - 0.0560?/ 2 + 0.0280a; + 0.2194?/ - 6.336 

x 2 - 0.4265?/ 3 + 6.118?/ 2 - 14.559a; - 32.047?/ + 77.711 

xy 2 + 0.0114y 3 - 0.1372?/ 2 + 0.0686x- 5.463y- 0.8231 

y 4 - UA77y 3 + 76. 724y 2 - 14.862a; - 188.419y + 214.345 

a;y 3 + 0.0280y 3 - 6.336y 2 + 0.1680a; + 1.316?/ - 2.016 



6 Appendix 

In this section we present the proof of Proposition 14.21 Proof: Let n E N; first 
we prove a result on the homogeneous components of degrees and 1 of the 
inverse of a square matrix A S Mat„ i „(i 7 '). 

Let A be a non singular matrix in M.&t n>n (F) and let B be the inverse of A. 
The homogeneous components Bq and B\ of B are given by 

B = A^ 1 Bt - -A Mx^ 1 (7) 

Define A A = J2i>2 ^ anc ^ ~ Si>2 -^«> so we nave 

A = A Q + A 1 + AA and B = B Q + B x + AB 
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Since AB = I, where I is the n x n identity matrix, we have 

(A + A X + AA)(B Q +B 1 + AB)=I 
and our claim is immediate. 

Now we prove the result of the proposition. Applying the classical least squares 
method to the linear system Mx ««we obtain 

x = (M t M)~ l M t v (8) 
p = v-Mx (9) 

Applying to © the homogeneous degree decomposition up to degree 1 we have 

PO + Pi = (VQ - MqXq) + (Vi - MqXx - MxXq) 

thus |3]) follows. 

Since (M*M) = M*M and (M t M) 1 = M*Mi + M(M , from formula © we 
have to first order, 

(M*M) _1 (M^Mo)- 1 - (MqMq) ~ 1 (MqMi + M t 1 M ){M^M )~ 1 

Up to degree 1, formula ([8]) becomes 

x + x x = (M* : M)o 1 (M^ + Mo'ui + Affuo) - (M* Af)^ 1 M*v 

= {M^Moy 1 (Mfro + M&t + MIvq - {MlM t + Af 1 t M )(M t M )" 1 M*w ) 

and so 

x = {MlM^Mlw 

Xl = (M^Mo)" 1 (Mfci + M\v - M^MxXa - M{M x ) 
thus the proof is concluded. 
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